Half-integer anomalous currents in 2D materials from a QFT viewpoint

Charge carriers in Dirac/Weyl semi-metals exhibit a relativistic-like behavior. In this work we propose a novel type of intrinsic half-integer Quantum Hall effect in 2D materials, thereby also offering a topological protection mechanism for the current. Its existence is rooted in the 2D parity anomaly, without any need for a perpendicular magnetic field. We conjecture that it may occur in disturbed honeycomb lattices where both spin degeneracy and time reversal symmetry are broken. These configurations harbor two distinct gap-opening mechanisms that, when occurring simultaneously, drive slightly different gaps in each valley, causing a net anomalous conductivity when the chemical potential is tuned to be between the distinct gaps. Some examples of promising material setups that fulfill the prerequisites of our proposal are also listed to motivate looking for the effect at the numerical and experimental level.

www.nature.com/scientificreports/ equalling e 2 2 . The A-QHE has triggered an intense research activity [13][14][15] , since it would be intrinsic to the material without the need for a (strong) external magnetic field as for the normal QHE, and when a 2D material with A-QHE is put in heterostructure with a superconducting material, the ensuing proximity effect could lead to the formation of the elusive Majorana fermions 16,17 , with great potential for (topologically protected against decoherence) qubits and thence quantum computation therewith. Interestingly enough, the inspiration behind also comes from a quantum field theory (QFT) viewpoint, see 18 . The intrinsic nature of the A-QHE is important in the latter context, as strong magnetic fields required for the "normal" QHE (which is non-intrinsic), would be destructive for the superconductor.

Theoretical setup from a QFT viewpoint
Setting the stage. We explore a honeycomb lattice with broken sublattice symmetry, as represented in Fig. 1b. It has been shown that the difference in the energy of the electrons belonging to different sublattices can be parameterized at the Hamiltonian level via an effective standard fermion mass for two two-component spinors, ±m , with opposite sign for each valley 19,20 . This induces a gap opening that preserves time reversal symmetry. In this case both valleys present symmetric gaps as represented in Fig. 1c. It is also known that the anomalous T-odd piece of the electric current generated by an external field acting on this system is related to a topological Chern-Simons term 21 . In other words, it depends solely on the sign of the mass term and not on its magnitude. In presence of an external gauge field in two spatial dimensions, one gets indicating that the net current vanishes when one sums up the contribution of both valleys 19 .
In this work we add two new ingredients to the case previously considered. First of all, we allow for arbitrary values of the gaps associated to each Dirac point, m k and m k ′ , enabling configurations where m k = m k ′ . This corresponds to a band structure where the two valleys have asymmetric gaps, as represented in Fig. 1d. Secondly, we include a chemical potential in order to explore certain regions of the band structure. The reason for this flexibility and how we calibrate the chemical potential will soon become clear.
For pure honeycomb lattices composed by one type of atom per lattice site, such as pristine graphene, the Lagrangian obtained from the continuum limit of the tight-binding model is very similar (apart from the Fermi velocity breaking Lorentz invariance) to the fermion sector of Quantum Electrodynamics in two space and one time dimensions 8, 9 , Here the index s labels the (electron) spin and the covariant derivative D α = ∂ α − (ie/ c)A α . A α is the gauge field associated to the electromagnetic interaction and the fields ψ s ≡ ψ s (t, r) denote four-component spinors that account for both valley and pseudospin index.
We choose to work in the Weyl basis,  www.nature.com/scientificreports/ To make the valley degree of freedom more transparent we rewrite our Lagrangian using the pseudo-chiral projections. Decomposing the fermion field as we can split the Lagrangian in Eq. (2) in two terms, associated to left and right pseudo-chiralities. We introduce untied gaps for each pseudo-chirality, allowing each projection to have different masses. The resulting Lagrangian is 22,23 where we included the chemical potential, µ , and a more general mass gap structure with m 2 s,k = m 2 s,k ′ . Without loss of generality, we may assume m k m ′ k > 0.
The fermion mass sector in two-component spinor language. For later usage, it is useful to rewrite the Lagrangian (5) a bit. The mass terms are written in four-component spinor language in the pseudo-chiral basis, meaning that (we suppress the spin index s for now.) using two two-component spinors ψ ± . The corresponding mass sector can be rewritten as keeping in mind that in (2 + 1)-dimensions, σ z plays the role of γ 0 in an appropriate basis 24 . This shows that we actually have two two-components spinors with opposite sign standard masses. This plays a pivotal role in our analysis.
Linear response theory applied to reduced quantum electrodynamics. Linear response theory investigates the reaction of a system under a small external influence assuming that this reaction can be studied at linear order in the external stimulus. Here we apply this procedure to calculate the conductivity of the just introduced continuum quantum field theory description of a material sample organized in a (2 + 1)-dimensional honeycomb lattice when a constant external electric field is applied in-plane, let's say � E = E� e y . The electric current we are interested in is then given by also in-plane but perpendicular to the applied E.
A certain care is needed since the underlying system is actually a mixed-dimensional theory: the charge carriers are constrained to a plane while the gauge fields live in a bulk. The theory to describe such systems has been developed previously and is known as pseudo-Quantum Electrodynamics (PQED) 25,26 or reduced-Quantum Electrodynamics (RQED) 27 . Here we present in detail how to introduce the interaction with the external electric field. For the sake of clarity, we refer in this subsection to the dimension of the system using the usual quantum field theory notation (n-space + 1-time)-dimension, in contradistinction with the condensed matter conventional notation used elsewhere in the manuscript.
The Kubo formalism has been successfully applied to several transport phenomena in Dirac/Weyl materials 28,29 . To extract the DC conductivity, we will benefit of some of the special properties of our mixed-dimensional system. It is a common trick to temporarily consider a time-varying electric field, and ultimately the DC limit ω → 0 will yield the desired constant field. The external (classical) electric field, and conveniently chosen vector potential, read in natural units where = c = 1 . Notice that this field configuration does not solve the classical Maxwell equations and cannot be considered a physical disturbance, unless for ω → 0 . We can introduce, however, a magnetic field in order to solve the Maxwell equations, but it will vanish in the limit ω → 0 , so we can omit it.
Having introduced the classical background, we proceed with the action in (3 + 1)-dimensions, given by since the (fermion) electronic degrees of freedom propagate on the z = 0 surface.
One can check the full gauge invariance of the (3 + 1)-dimensional action in Eq. (11), applying with 4-vector x = (� x, z) . We tacitly assume so-called absolute boundary conditions 30 0,x,y z=0 = j 0,x,y , as also used in 31,32 . Notice that these are consistent with using the Landau gauge. The Noether current is easily identified as This current is gauge invariant and conserved on-shell, as expected. It is also independent from the z-coordinate.
The third component j z can be written down formally, but it does not correspond to the transport of any physical electric charge, as can be easily checked using the integrated Gauss' law. So effectively, j z = 0. We notice that j µ has mass dimension 2, as corresponding to a (2 + 1)-dimensional current. Let us introduce its (3 + 1)-dimensional version, to properly disturb the action S with to set up the appropriate Kubo linear response theory.
The (still mixed-dimensional) perturbed action can thus be written as Following e.g. 32 and transforming to momentum (Fourier) space, it can be shown that the action in Eq. (14) can be reduced to an equivalent, be it purely (2 + 1)-dimensional, action We recognize here RQED, supplemented with an external disturbance. Doing so, we can now depart from this consistent (2 + 1)-dimensional action to extract the DC conductivity via a Kubo relation. There is no more reference to the obsolete z-direction. This also makes clear why we chose the vector potential (10), which is by no means unique, as this choice allows for a full dimensional reduction of the relevant dynamics. Closely following 33 , we have a disturbed Hamiltonian given by H = � j · � a , from now on always assuming z = 0 . The expectation value of the current to leading order in a becomes or concretely, after using time translational invariance, we get the Kubo relation for the DC anomalous conductivity, As we are interested in the ω → 0 limit, it will be sufficient to evaluate the integral in Eq. (17) up to O(ω) . As the integrand in the r.h.s. of Eq. (17) is actually the retarded correlator, we can use the fact that retarded and Euclidean self-energy correlator coincide at ω = 0 . We may use (Ref. 34 , Eq. (2.11) to find they do so up to O(ω 2 ) , an identification valid since the Euclidean photon self-energy ij has no pole at zero frequency (see Eq. (22)). In fact, from Eq. (20) it follows the spectral function will only be non-vanishing from a non-zero threshold onwards, exactly due to the infrared being protected by the mass scale.
The RQED self-energy at zero momentum. Assuming first a standard Dirac mass m for a single twocomponent spinor and a chemical potential µ , it can be shown that the low-energy limit of the induced T-odd part in the RQED photon self-energy is one-loop exact 35 . This is the Coleman-Hill theorem 36 generalized to RQED, resulting in a "topological" Chern-Simons mass term in the effective action. Moreover, that relevant contribution is also independent of the Fermi velocity v F 35 . At one-loop the (transverse and gauge invariant) polarization tensor is only depending on fermion propagators, which are just the same as in usual planar QED. www.nature.com/scientificreports/ The required (Euclidean) one-loop photon self-energy-which corresponds to the current-current correlator-reads where in general k = (k 0 + iµ, v F k x , v F k y ).
A non-vanishing contribution to ˜ xy can only arise from the T-odd piece of the above integral, a piece which is necessarily proportional to the Levi-Civita tensor ε ijk . The latter comes exclusively from the product of three γ-matrices, based on Tr(γ i γ j γ k ) = −2ε ijk . So we get, with p 0 ≡ ω, where Given the explicit p k in front of expression (19) and the eventual ω → 0 limit, it is sufficient for our purposes to compute π(0), the last step based on the residue theorem, the integral vanishes if m 2 < µ 2 as the contour can then be chosen to not encircle any of the poles occurring at p 0 = −iµ ± i|m| . Putting everything together, we actually have This result coincides with that of e.g. [37][38][39][40] , where the temperature T → 0 limit is understood whenever necessary. Let us refer the interested reader to e.g. 41,42 for alternative derivations for several types of Quantum Hall conductivities.
To avoid any doubts about the DC limit, we find it useful to explicit verify that our result is independent on how the zero momentum limit is taken, namely p 0 = 0 , lim p x,y → 0 vs. p x,y = 0 , lim p 0 → 0 . Adopting the Feynman trick, or, after substituting ℓ 0 = q 0 + αp 0 , ℓ x,y = v F (q x,y + αp x,y ), As the net dependence is clearly on � p 2 = p 2 0 + v 2 F p 2 x + v 2 F p 2 y , from the above expression it is already obvious the aforementioned limits do commute. For the momentum integration, we can follow exactly the same procedure as in the � p 2 = 0 case. So, after integrating over the Feynman parameter, we end up with It can be verified that in the µ → 0 limit, this result coincides with the one quoted in e.g. 24  for a single massive two-component spinor. This DC conductivity will receive no further corrections, given the one-loop exactness of � yx (ω) for vanishing frequency. This is suggestive of the fact the conductivity might have a topological origin, something that will be discussed further in the next subsection. To derive our final expression for the anomalous conductivity, in the case that only a single spin-band contributes, albeit two different two-component spinors (one per Dirac point), www.nature.com/scientificreports/ Assuming now that m k > 0 , m k ′ > 0 and m 2 k ′ < µ 2 < m 2 k , the predicted value for the conductivity becomes which corresponds to a half-integer anomalous quantum Hall conductivity capable of sustaining a current along an externally applied electric field. A pivotal role is played here by assuming the chemical potential to be in the gap between the two gaps, otherwise we would get a net cancellation of the two terms in Eq. (27) yet again.
Reintroducing the spin degree of freedom, we have where m s,k and m s,k ′ are now the gaps associated to each valley, per spin choice. This is our main result. We note that for the configuration considered in previous works 12, 19 , where |m s,k | = |m s,k ′ | , the net current vanishes as expected, but let us consider the case where the gaps differ, represented in Fig. 1d. Then the system does enjoy neither a space inversion nor time reversal invariance, implying an uplifting of the pseudo-chiral fermions degeneracy. The most immediate way to achieve such a configuration is to work simultaneously with a double gap opening mechanism. On top of the inversion symmetry breaking mass term, δτ z ⊗ σ z -with τ designating the valley whilst σ relating to sublattice-discussed previously 19 , we propose to add an interaction capable of inducing a time reversal symmetry breaking. This corresponds to a mass term where the sign of the masses remain the same in both valleys, �τ 0 ⊗ σ z . A configuration subject to both symmetry breaking mechanisms will lead to different gaps in each valley as required.
The coupling of the masses to spin must be chosen carefully, otherwise currents associated to different spins cancel each other. First of all, we require that a spin flip causes a global change of sign, not a relative one. This is crucial otherwise for each spin a different valley will be favored, inflicting a cancellation. That is the reason a Kane-Mele type of configuration 43 does not satisfy our conditions. In order to avoid this situation we choose to add the following mass couplings to the Lagrangian, with � > δ.
Notice however that a flip in spin induces a flip in mass sign. We can check from Eq. (29) that this also would induce a cancellation between currents belonging to the same valley with different spin. We will consider this observation below, when we seek for suitable materials that fulfill our conditions. Synthesizing the discussion above, a planar material capable to support our proposed half-integer A-QHE must (i) contain Dirac points (ii) present inversion symmetry breaking (iii) present time reversal symmetry breaking and (iv) present a split in spin bands.
If these conditions are met, we thus predict from Eq. (29) a DC current given by where we reintroduced the units.
Compared to the Haldane model 13 , there are two differences. We consider two symmetry breakings-apart from a spin polarization to lift the spin degeneracy, which is also present in 13 -namely time reversal and space inversion symmetry, while in the Haldane model only time reversal is broken whilst space inversion invariance is preserved. Different from 13 , we do not propose band inversion but different gaps per valley, which is a crucial feature to allow for a chemical potential in the gap between the gaps.

Half-integer topology behind the anomalous conductivity
As we explained before, the anomalous conductivity (30) will receive no further quantum corrections, which is suggestive of a topological protection mechanism. That this is indeed the case will be scrutinized now. We lay the connection with a half-integer A-QHE in detail, not unlike the situation of the magnetically doped (BiSb) 2 (TeSe) 2 materials reviewed in 14 albeit that the underlying mechanism is quite different in our case. Indeed, we are not looking at the single Dirac cone per surface of a three-dimensional topological insulator, rather our judicious choice of model parameters, our genuine two-dimensional material eliminates one of the two cones from the game. Moreover, we do not have to rely on band inversion. What is common in both situations is that the effective low energy dynamics around the relevant Dirac point is attributable to a continuum Chern-Simons term, radiatively generated by the massive fermions.
For the benefit of the reader, let us briefly recapitulate the reasoning, see e.g. 33 . Reconsidering Eq. (17), and using the undisturbed energy eigenbasis, |n� , and Hamiltonian time evolution with initial time t = 0 , � j(t) = e i H 0 t � je −i H 0 t , we may rewrite it as (27) � j = σ xy E� e y , σ xy = 1 2 www.nature.com/scientificreports/ Integrating over τ and expanding to leading order in ω again, we arrive at No divergence in 1 ω hampers the discussion as the photon self-energy (viz. current correlator) starts at O(ω) , see 35 for an explicit proof based on a Ward identity (in se gauge invariance).
Proceeding as in 33 by introducing for both α filled bands and β unfilled bands the Bloch wave functions, � α,β � p (� x) = e i� p·� x u α,β � p (� x) , we can rewrite Eq. (32) using a compact notation where the 2-dimensional p lives on the Brillouin zone T 2 . This expression can be further massaged into once the current in Fourier space is identified as � j = e ∇ � p H , with H = e −i� p·� x He i� p·� x , the Hamiltonian with momentum shifted over p . To arrive at Eq. (34), it was tacitly assumed that the bands are gap-separated, with chemical potential µ in such gap. The final outcome is the well-known sum over the integer Chern numbers C α per filled band 28 , thereby disposing the topological character of the conductivity.
However, we should be careful when applying Eq. (34) directly to our case, as one of the underlying conditions is not met. Indeed, a crux of our setup is that the chemical potential µ was chosen in the gap between the 2 Dirac points. Evidently, this leads to a "cutting" of the respective band, see Fig. 1d. Let us therefore analyze the expression (36) a bit better in the case the underlying low energy effective Hamiltonian per valley per spin is of the type 44 The few materials that we will identify in the next section, with the desired gap structure, are all of this type. To be more precise, they are of the graphene-type, but supplemented with a specific mass term along the σ z -direction.
Then, using the "ray method", it was shown in 44 that each massive Dirac point (cone) contributes with ± 1 2 "topological charge". If the full (smooth) band is taken into account for a double cone system (as we have), this will automatically lead to an integer Chern number and thence σ = ± e 2 h , 0 as possible conductivities. However, due to our judicious choice of µ in the finite intervalley gap between the K and K ′ Dirac points, we are effectively cutting away part of that band and only a single cone will contribute to the conductivity. Denoting with α ′ all completely filled bands and singling out that "partially cut" band (the encircled piece), we actually have where the encircled piece is only there for the spin projection that is partially filled, the sign in front depending on which of the 2 Dirac points is the contributing one. Remember also that we lifted the spin degeneracy to avoid a net cancellation of the current; the other spin projection band is thus not filled and as such not contributing. This encircled piece exactly confirms the result (28) found before using an explicit quantum field theory computation and as promised, despite its half-integer character, the above reasoning shows it still carries a topological meaning. Indeed, small perturbations of the band geometry will preserve the fact that only a single Dirac point contributes to the Berry curvature integral, which always happens with a value of ± 1 2 . This being said, the cutting of the band is a non-smooth operation, leading to a non-smooth integration zone, explaining why strictly speaking, we dot not find a topological integer.
In future work, we will analyze the edge mode spectrum in full detail by solving the gapped Dirac equation associated to Eq. (5) for zero modes on a semi-infinite strip, but the judicious choice of the chemical potential in the gap between the different gaps, is expected to also lead to a halving of the number of (normalizable) edge modes, so the bulk-edge correspondence 45 , stills holds, thereby offering an alternative view on the origin of the half-integer conductivity and its robustness. Although in a different context-and using different conventions, and not necessarily looking at zero modes-Ref. 46 , Eq. (5) already illustrates that the presence of a chemical potential can lead to one normalizable and one non-normalizable mode.

Promising materials
A few candidates present the features that meet the criteria laid out in Sect. II, as calculated using ab initio methods. In general, monolayer crystals with high potential to be applied in valleytronics physics are good candidates. First of all, manganese chalcogenophosphate MnPSe 3 is organized in honeycomb lattices and it is intrinsically gapped. The material counts on two mechanisms to open the gaps: the antiferromagnetic order coupled to the sublattice and a spin-valley coupling, being responsible for time symmetry and inversion breaking, respectively 47,48 . The situation corresponds exactly to the masses represented in Eq. (30). However, the spin remains degenerate due to the antiferromagnetic (AFM) coupling between two inequivalent Mn ions. www.nature.com/scientificreports/ In order to generate the necessary split between the spin bands, two types of setups can be prepared. One of them consists of doping the material, replacing for example one Mn atom by a Zn atom. The imbalance in spin will generate a magnetization, and consequently a spin splitting in the bands, while still retaining the different valleys 49 , as can be seen in Fig. 2. Note that we have the same structure schematically represented in Fig. 1e with a spin up (down) valley splitting of 33 (40) meV and a spin splitting of 89 (93) meV at the K' (K) point. The band structure in Fig. 2 was obtained via a density functional theory (DFT) calculation with the VASP package 50 for the same system represented in Ref. 49 , Fig. 2a. Exchange and correlation are described within the Perdew-Burke-Ernzerhof functional 51 with a on site Coulomb interaction of 5 eV applied to the Mn atoms ( GGA + U) 52 .
Another possible way is to use an heterojunction of different 2D materials such as MnPSe 3 /CrBr 3 53 , MnPSe 3 /MoS 2 54 and WS 2 /h-VN 55 . Besides the aforementioned heterostructures, dichalcogenides like NbSe 2 and WS 2 also seem promising to simultaneously lift both valley and spin degeneracy. This can be either reached by using a substrate 56 to WS 2 , or as by using a 2D magnetic semiconductor as monolayer NbSe 2 with its large valley-polarized state 57 . In such system, no external influence by either doping, substrating or van der Waals couplings is necessary, rather the magnetic state is sufficient to completely split the bands, resulting in the generation of a valley-polarized state with a spontaneously occurring valley current 58,59 . Our proposed effect thus is a specific type of a spin-polarized Quantum valley Hall effect, where the novel half-integer nature is due to a fine-tuning of the chemical potential. The latter can, for example, be realized in practice by applying a suitable external voltage.

Conclusion
We have proposed a half-integer version of the A-QHE, intrinsic to certain specifically gapped 2D materials, with as raison d'être the parity anomaly.
We offered several setups that follow the criteria to source this mechanism. Therefore, we expect that our findings may open a new window to explore anomaly driven current in planar materials, offering potential advantages in the search for new anomalous transport phenomena.
Interestingly, for another class of transition metal chalcogenides, like V 2 O 3 , recent computations of 15,60 suggest that the two valleys contribute with equal sign mass to the dynamical Chern-Simons term, something which should be related to the strong spin-orbit coupling. This could lead to an intrinsic integer rather than half-integer conductivity. This deserves further scrutiny in the future.